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The range of expected amplitudes and spectral slopes of relic (squeezed) 
gravitational waves, predicted by theory and partially supported by obser- 
vations, is within the reach of sensitive gravity-wave detectors. In the most 
favorable case, the detection of relic gravitational waves can be achieved by 
the cross-correlation of outputs of the initial laser interferometers in LIGO, 
VIRGO, GEO600. In the more realistic case, the sensitivity of advanced 
ground-based and space-based laser interferometers will be needed. The spe- 
cific statistical signature of relic gravitational waves, associated with the phe- 
nomenon of squeezing, is a potential reserve for further improvement of the 
signal to noise ratio. 



I. INTRODUCTION 

It is appropriate and timely to discuss the detection of relic gravitatational waves at 
the experimental meeting like this one. We are in the situation when the advanced laser 
interferometers, currently under construction or in a design phase, can make the dream of 
detecting relic gravitons a reality. The detection of relic gravitational waves is the only way 
to learn about the evolution of the very early Universe, up to the limits of Planck era and 
Big Bang. 

The existence of relic gravitational waves is a consequence of quite general assumptions. 
Essentially, we rely only on the validity of general relativity and basic principles of quantum 
field theory. The strong variable gravitational field of the early Universe amplifies the in- 
evitable zero-point quantum oscillations of the gravitational waves and produces a stochastic 
background of relic gravitational waves measurable today []]]. It is important to appreci- 
ate the fundamental and unavoidable nature of this mechanism. Other physical processes 
can also generate stochastic backgrounds of gravitational waves. But those processes ei- 
ther involve many additional hypotheses, which may turn out to be not true, or produce a 
gravitational wave background (like the one from binary stars in the Galaxy) which should 
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be treated as an unwanted noise rather than a useful and interesting signal. The scientific 
importance of detecting relic gravitational waves has been stressed on several occasions (see, 
for example, 0- ||]). 

The central notion in the theory of relic gravitons is the phenomenon of superadiabatic 
(parametric) amplification. The roots of this phenomenon are known in classical physics, 
and we will remind its basic features. As every wave-like process, gravitational waves use the 
concept of a harmonic oscillator. The fundamental equation for a free harmonic oscillator is 

q + u 2 q = 0, (1) 

where q can be a displacement of a mechanical pendulum or a time- dependent amplitude 
of a mode of the physical field. The energy of the oscillator can be changed by an acting 
force or, alternatively, by a parametric influence, that is, when a parameter of the oscillator, 
for instance the length of a pendulum, is being changed. In the first case, the fundamental 
equation takes the form 

q + U 2 q = f(t), (2) 
whereas in the second case Eq. ([]]) takes the form 

q + LU 2 (t)q = 0. (3) 

Equations (0) and (|3|) are profoundly different, both, mathematically and physically. 

Let us concentrate on the parametric influence. We consider a pendulum of length L 
oscillating in a constant gravitational field g. The unperturbed pendulum oscillates with the 
constant frequency to = J~gfL. Fig. la illustrates the variation of the length of the pendulum 
L(t) by an external agent, shown by alternating arrows. Since L(t) varies, the frequency 
of the oscillator does also vary: uj(t) = \Jg/L(t). The variation L(t) does not need to be 
periodic, but cannot be too much slow (adiabatic) if the result of the process is going to be 
significant. Otherwise, in the adiabatic regime of slow variations, the energy of the oscillator 
E and its frequency uj do change slowly, but E/u remains constant, so one can say that the 
"number of quanta" E/Tioj in the oscillator remains fixed. In other words, for the creation of 
new "particles - excitations" , the characteristic time of the variation should be comparable 
with the period of the oscillator and the adiabatic behaviour should be violated. After some 
duration of the appropriate parametric influence, the pendulum will oscillate at the original 
frequency, but will have a significantly larger, than before, amplitude and energy. This is 
shown in Fig.lb. Obviously, the energy of the oscillator has been increased at the expense 
of the external agent (pump field). For simplicity, we have considered a familiar case, when 
the length of the pendulum varies, while the gravitational acceleration g remains constant. 
Variation of g would represent a gravitational parametric influence and would be even in a 
closer analogy with what we study below. 

A classical oscillator must have a non-zero initial amplitude for the amplification mech- 
anism to work. Otherwise, if the initial amplitude is zero, the final amplitude will also be 
zero. Indeed, imagine the pendulum strictly at rest, hanging stright down. Whatever the 
variation of its length is, it will not make the pendulum to oscillate and gain energy. In 
contrast, a quantum oscillator does not need to be excited from the very beginning. The 
oscillator can be initially in its vacuum quantum-mechanical state. The inevitable zero-point 
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quantum oscillations are associated with the vacuum state energy ~%lo. One can imagine a 
pendulum hanging stright down, but fluctuating with a tiny amplitude determined by the 
"half of the quantum in the mode". In the classical picture, it is this tiny amplitude of 
quantum-mechanical origin that is being parametrically amplified. 




Fig. 1. Parametric amplification, a) variation of the 
length of the pendulum, b) increased amplitude of 
oscillations. 

The Schrodinger evolution of a quantum oscillator depends crucially on whether the 
oscillator is being excited parametrically or by a force. Consider the phase diagram (q,p), 
where q is the displacement and p is the conjugate momentum. The vacuum state is described 
by the circle in the center (see Fig. 2). The mean values of q and p are zeros, but their 
variances (zero-point quantum fluctuations) are not zeros and are equal to each other. Their 
numerical values are represented by the circle in the center. Under the action of a force, 
the vacuum state evolves into a coherent state. The mean values of p and q have increased, 
but the variances are still equal and are described by the circle of the same size as for the 
vacuum state. On the other hand, under a parametric influence, the vacuum state evolves 
into a squeezed vacuum state. [For a recent review of squeezed states see, for example, 
H and references there.] Its variances for the conjugate variables q and p are significantly 
unequal and are described by an ellipse. As a function of time, the ellipse rotates with 
respect to the origin of the (q, p) diagram, and the numerical values of the variances oscillate 
too. The mean numbers of quanta in the two states, one of which is coherent and another 
is squeezed vacuum, can be equal (similar to the coherent and squeezed states shown in 
Fig. 2) but the statistical properties of these states are significantly different. Among other 
things, the variance of the phase of the oscillator in a squeezed vacuum state is very small 
(squeezed). Graphically, this is reflected in the fact that the ellipse is very thin, so that that 
the uncertainty in the angle between the horizontal axis and the orientation of the ellipse is 
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very small. This highly elongated ellipse can be regarded as a portarait of the gravitational 
wave quantum state that is being inevitably generated by parametric amplification, and 
which we will be dealing with below. 



coherent state 




Fig. 2. Some quantum states of a harmonic oscillator. 

A wave-field is not a single oscillator, it depends on spatial coordinates and time, and may 
have several independent components (polarization states). However, the field can be de- 
composed into a set of spatial Fourier harmonics. In this way we represent the gravitational 
wave field as a collection of many modes, many oscillators. Because of the nonlinear charac- 
ter of the Einstein equations, each of these oscillators is coupled to the variable gravitational 
field of the surrounding Universe. For sufficiently short gravitational waves of experimental 
interest, this coupling was especially effective in the early Universe, when the condition of 
the adiabatic behaviour of the oscillator was violated. It is this homogeneous and isotropic 
gravitational field of all the matter in the early Universe that played the role of the external 
agent - pump field. The variable pump field acts parametrically on the gravity-wave oscilla- 
tors and drives them into multiparticle states. Concretely, the initial vacuum state of each 
pair of waves with oppositely directed momenta evolves into a highly correlated state known 
as the two- mode squeezed vacuum state ||, 0. The strength and duration of the effective 
coupling depends on the oscillator's frequency. They all start in the vacuum state but get 
excited to various amounts. As a result, a broad spectrum of relic gravitational waves is 
being formed. This spectrum is accessible to our observations today. 

Let us formulate the problem in more detail. 

II. COSMOLOGICAL GRAVITATIONAL WAVES 

In the framework of general relativity, a homogeneous isotropic gravitational field is 
decribed by the line element 
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ds 2 = c 2 dt 2 - a 2 (t)5 ij dx i dx j = a 2 {7])[drf - %dxW]. 



(4) 



In cosmology, the function a(t) (or a{rf)) is called scale factor. In our discussion, it will 
represent gravitational pump field. 

Cosmological gravitational waves are small corrections to the metric tensor. They 
are defined by the expression 



ds 2 = a 2 (r})[dr} 2 - {5 tJ + h^dx'dx 3 }. 



(5) 



The functions hij(r), x) can be expanded over spatial Fourier harmonics e mx and e _mx , 
where n is a constant wave vector. In this way, we reduce the dynamical problem to the 
evolution of time-dependent amplitudes for each mode n. Among six functions there are 
only two independent (polarization) components. This decomposition can be made, both, 
for real and for quantized field hij. In the quantum version, the functions are treated 
as quantum-mechanical operators. We will use the Heisenberg picture in which the time 
evolution is carried out by the operators while the quantum state is fixed. This picture 
is fully equivalent to the Schrodinger picture, discussed in the Introduction, in which the 
vacuum state evolves into a squeezed vacuum state while the operators are time independent. 
The Heisenberg operator for the quantized real field can be written as 



C f°° 2 s 1 



hn{v)e in * c n + h n (v)e 



(6) 



where C is a constant which will be discussed later. The creation and annihilation operators 

satisfy the conditions [c n , c m ] = 5 s r s 5 3 (n — m), c n |0) = 0, where |0) (for each n and s) 
is the fixed initial vacuum state discussed below. The wave number n is related with the 
wave vector n by n = (SijU 1 ^) 1 ^ 2 . The two polarization tensors Pij(n) (s = 1,2) obey the 
conditions 



P ijn j = 0, Pij5 ij = 0, PijP ij = 25 ss ,, Pij(-n) = P i:j (n). 
The time evolution, one and the same for all n belonging to a given n, is represented by 

s 

the complex time-dependent function h n (v)- This evolution is dictated by the Einstein 

s 

equations. The nonlinear nature of the Einstein equations leads to the coupling of hn{V) 
with the pump field a(i]). For every wave number n and each polarization component s, the 

functions hn(f]) have the form 



hn{v) = ^[Un(v)+Z(v)}, 



(7) 



where u n (r]) and v n {rf) can be expressed in terms of the three real functions (the polarization 
index s is omitted): r n - squeeze parameter, <p n - squeeze angle, 6 n - rotation angle, 



u n = e 



cosh 



= e - i(e "- 2 ^ } sinh r„. 



(8) 



The dynamical equations for u n (r]) and v n (rj) 
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■ du n .a' „ .dv n .a' 

i—— = nu n + i—v, %—- = nv n + i— u n (9) 
drj a drj a 

lead to the dynamical equations governing the functions r n (rf), <fi n (v)-> 0n(ij) 01 : 

o! o! q! 

r' n = — cos20 n , (p' n = — n sin20 n coth 2r n , 9'——n sin20 n tanh r n , (10) 

a a a 

where ' = d/drj, and the evolution begins from r n = 0. This value of r n characterizes 
the initial vacuum state |0) which is defined long before the interaction with the pump field 
became effective, that is, long before the coupling term a'/ a became comparable with n. The 
constant C should be taken as C = Vl6vr Ipi where Ipi = (Gh/c 3 ) 1 ^ 2 is the Planck length. 
This particular value of the constant C guarantees the correct quantum normalization of 
the field: energy ^Tiuj per each mode in the initial vacuum state. The dynamical equations 
and their solutions are identical for both polarization components s. 

Equations @ can be translated into the more familiar form of the second-order differ- 

S g g * s 

ential equation for the function fi> n (r]) = u n {rj) + v n (r)) = a(r])h n (v) : 



2 a" 
n z 



0. (11) 



Clearly, this is the equation for a parametrically disturbed oscillator (compare with Eq. 
(H)). In absence of the gravitational parametric influence represented by the term a" /a, the 
frequency of the oscillator defined in terms of 77-time would be a constant: n. Whenever the 
term a" /a can be neglected, the general solution to Eq. (|TT|) has the usual oscillatory form 

fi n ( V ) = A n e~ m i + B n e m \ (12) 

where the constants A n , B n are determined by the initial conditions. On the other hand, 
whenever the term a" /a is dominant, the general solution to Eq. ([11]) has the form 

~aJ' ^ 

In fact, this approximate solution is valid as long as n is small in comparison with \a'/a\. 
This is more clearly seen from the equivalent form of Eq. (0) written in terms of the 
function h n {rf) |5|: 

ti^ + 2^ti n + n 2 h n = 0. (14) 

For growing functions a(r]), that is, in expanding universes, the second term in Eq.(|13|) is 
usually smaller than the first one (see below), so that, as long as n C a' /a, the dominant 
solution is the growing function fi n (v) = C n a(rj), and 

h n = const. (15) 

Equation ([□]) can be also looked at as a kind of the Schrodinger equation for a particle 
moving in presence of the effective potential U(rj) = a" /a. In the situations that are normally 
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considered, the potential U{rj) has a bell-like shape and forms a barrier (see Fig. 3). When 
a given mode n is outside the barrier, its amplitude h n is adiabatically decreasing with 
time: h n oc e , > . This is shown in Fig.3 by oscillating lines with decreasing amplitudes of 
oscillations. The modes with sufficiently high frequencies do not interact with the barrier, 
they stay above the barrier. Their amplitudes h n behave adiabatically all the time. For 
these high-frequency modes, the initial vacuum state (in the Schrodinger picture) remains 
the vacuum forever. On the other hand, the modes that interact with the barrier are subject 
to the superadiabatic amplification. Under the barrier and as long as n < a! j a, the function 
h n stays constant instead of the adiabatic decrease. For these modes, the initial vacuum 
state evolves into a squeezed vacuum state. 




Fig. 3. Effective potential U(r|). 

After having formulated the initial conditions, the present day behaviour of r n , <f) n , 6 n 
(or, equivalently, the present day behaviour of h n ) is essentially all we need to find. The 
mean number of particles in a two-mode squeezed state is 2 sinh 2 r n for each s. This number 
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determines the mean square amplitude of the gravitational wave field. The time behaviour 
of the squeeze angle (fi n determines the time dependence of the correlation functions of the 
field. The amplification (that is, the growth of r n ) governed by Eq. ( |i0|) is different for 
different wave numbers n. Therefore, the present day results depend on the present day 
frequency v {y — cnjlixa) measured in Hz. 

In cosmology, the function H = a/a = ca'/a 2 is the time-dependent Hubble parameter. 
The function I = c/H is the time-dependent Hubble radius. The time-dependent wavelength 
of the mode n is A = lixa/n. The wavelength A has this universal definition in all regimes. In 
contrast, the v defined as v = cn/lixa has the usual meaning of a frequency of an oscillating 
process only in the short-wavelength (high-frequency) regime of the mode n, that is, in the 
regime where A /. As we have seen above, the qualitative behaviour of solutions to Eqs. 
'PI), ( |T4| ) depends crucially on the comparative values of n and a' /a, or, in other words, on 



the comparative values of A (77) and l(n). This relationship is also crucial for solutions to Eq. 
(PH|) as we shall see now. 

In the short-wavelength regime, that is, during intervals of time when the wavelength 
X(i]) is shorter than the Hubble radius l(r]) = a 2 /a', the term n in ( |T0| ) is dominant. The 
functions (fi n (v) an d 9n(v) are 0« = ~ n (v + Vn), @n = (fin where r\ n is a constant. The factor 
cos20 n is a quickly oscillating function of time, so the squeeze parameter r n stays practically 
constant. This is the adiabatic regime for a given mode. 

In the opposite, long-wavelength regime, the term n can be neglected. The function (fi n 
is tan (^(77) ~ const/ a 2 (n), and the squeeze angle quickly approaches one of the two values: 
(fin = or (fi n = ii (analog of "phase bifurcation" ||). The squeeze parameter r n (n) grows 
with time according to 

r n (v) ~ , (16) 

a* 

where a* is the value of a(n) at 77*, when the long- wavelength regime, for a given n, begins. 
The final amount of r n is 



r n w In— , (17) 
a* 

where a** is the value of a(n) at 77**, when the long-wavelength regime and amplification 
come to the end. It is important to emphasize that it is not a "sudden transition" from one 
cosmological era to another that is responsible for amplification, but the entire interval of 
the long-wavelength (non-adiabatic) regime. 

After the end of amplification, the accumulated (and typically large) squeeze parameter 
r n stays approximately constant. The mode is again in the adiabatic regime. In course of 

S S * 

the evolution, the complex functions u n (n) + v n {rf) become practically real, and one has 

S S * 

h n (v) ~ hniv) ~ c rn cos (fi n (jf) . Every amplified mode n of the field (gj) takes the form 
of a product of a function of time and a (random, operator-valued) function of spatial 
coordinates; the mode acquires a standing-wave pattern. The periodic dependence cos (fi n {n) 
will be further discussed below. 

It is clearly seen from the fundamental equations (|TUD, QTTJ) , fll^p that the final results 
depend only on a(rj). Equations do not ask us the names of our favorite cosmological 
prejudices, they ask us about the pump field a(r]). Conversely, from the measured relic 
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gravitational waves, we can deduce the behaviour of a(rj), which is essentially the purpose 
of detecting the relic gravitons. 



III. COSMOLOGICAL PUMP FIELD 

With the chosen initial conditions, the final numerical results for relic gravitational waves 
depend on the concrete behaviour of the pump field represented by the cosmological scale 
factor a{rj). We know a great deal about a(rj). We know that a(rj) behaves as a{rf) oc rf 
at the present matter-dominated stage. We know that this stage was preceeded by the 
radiation-dominated stage a(rj) oc 77. At these two stages of evolution the functions a{rf) are 
simple power-law functions of 77. What we do not know is the function a{rf) describing the 
initial stage of expansion of the very early Universe, that is, before the era of primordial 
nucleosynthesis. It is convenient to parameterize a {rf) at this initial stage also by power-law 
functions of 77. First, this is a sufficiently broad class of functions, which, in addition, allows 
us to find exact solutions to our fundamental equations. Second, it is known |TJ that the 
pump fields a(rf) which have power-law dependence in terms of 77, produce gravitational 
waves with simple power-law spectra in terms of v. These spectra are easy to analyze and 
discuss in the context of detection. 

We model cosmological expansion by several successive eras. Concretely, we take a{rf) 
at the initial stage of expansion (i-stage) as 

a(r ] )=l \v\ 1+ ' 3 , (18) 

where 77 grows from —00, and 1 + (5 < 0. We will show later how the available observational 
data constrain the parameters l and f3. The i-stage lasts up to a certain rj = 771, 771 < 0. To 
make our analysis more general, we assume that the i-stage was followed by some interval 
of the z-stage (z from Zeldovich). It is known that an interval of evolution governed by the 
most "stiff" matter (effective equation of state p — e) advocated by Zeldovich, leads to a 
relative increase of gravitational wave amplitudes |]J . It is also known that the requirement 
of conistency of the graviton production with the observational restrictions does not allow 
the "stiff" matter interval to be too much long ||1Q|| . However, we want to investigate any 
interval of cosmological evolution that can be consistently included. In fact, the z-stage of 
expansion that we include is quite general. It can be governed by a "stiffer than radiation" 



nj matter, as well as by a "softer than radiation" matter. It can also be simply a part 
of the radiation-dominated era. Concretely, we take 0(77) at the interval of time from rjx to 
some 1] S (z-stage) in the form 

a(77) = l a z {r l -r lp ) 1+p % (19) 

where 1 + (3 S > 0. For the particular choice j3 s = 0, the z-stage reduces to an interval of 
expansion governed by the radiation-dominated matter. Starting from rj s and up to 772 the 
Universe was governed by the radiation-dominated matter (e-stage). So, at this interval of 
evolution, we take the scale factor in the form 

a(rj) = l a e (rj - rj e ). (20) 

And, finally, from 77 = 772 the expansion went over into the matter-dominated era (m-stage): 
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a{rj) = l a m {v - r) m ) . (21) 

A link between the arbitrary constants participating in Eqs. fll8|) - (21) is provided by the 
conditions of continuous joining of the functions a(rj) and a'ijf) at points of transitions 771 , 

We denote the present time by r]f> (R from reception). This time is defined by the 
observationally known value of the present-day Hubble parameter H(j]r) and Hubble radius 
Ijj = c/H(j]r). For numerical estimates we will be using Ijj « 2 x 10 28 cm. It is convenient 
to choose fjR — r] m = 1, so that a(i] R ) = 2l H . The ratio 

aiVn)/^) = C2 
is believed to be around (2 = 10 4 . We also denote 



a(m)/a(Vs) = Cs , a{r) s ) / a^x) = Ci • 

With these definitions, all the constants participating in Eqs. ( |TB| ) - ( |2T| ) (except parameters 
/3 and /3 S which should be chosen from other considerations) are being expressed in terms of 
Ih, (2, Cs, and (±. For example, 

I I 11 + 11 

M = — — • 

2C2 2 C s Ci^ 

The important constant l Q is expressed as 

lo = bl H (^ 1 (!!(M, (22) 

where b = 2 2+/3 /|l + (3\ l+ ^ ' . Note that 6 = 1 for (3 = —2. [This expression for l Q may help to 
relate formulas written here with the equivalent treatment [ff^] which was given in slightly 
different notations.] The sketch of the entire evolution 0(77) is given in Fig A. 




Fig. 4. Scale factor a(r\). 
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We work with the spatially-flat models (Q). At every instant of time, the energy density 
e(r)) of matter driving the evolution is related with the Hubble radius l(rj) by 



Ke{r)) 



(23) 



where k = 87rG/c 4 . For the case of power-law scale factors a(r]) oc r/ +/3 , the effective pressure 
p(i]) of the matter is related with the e{rf) by the effective equation of state 



V 



1-/3 
3(1 + /?)' 



For instance, p - 
of the evolution 



Ofor/3 

T 



l,p=\e for /3 = 0, 



P 



-e for /3 



(24) 



-2, and so on. Each interval 



(pll) is governed by one of these equations of state. 
In principle, the function a (77) could be even more complicated than the one that we 
consider. It could even include an interval of the early contraction, instead of expansion, 
leading to the "bounce" of the scale factor. In case of a decreasing a{rj) the gravitational- wave 
equation can still be analyzed and the amplification is still effective . However, the Einstein 
equations for spatially-flat models do not permit a regular "bounce" of 0(77) (unless e vanishes 
at the moment of "bounce"). Possibly, a "bounce" solution can be realized in alternative 
theories, such, for example, as string-motivated cosmologies [IB 
spectral slopes of gravitational waves produced in "bounce" cosmologies, see WA . 



For a recent discussion of 



IV. SOLVING GRAVITATIONAL WAVE EQUATIONS 

The evolution of the scale factor a{j]) given by Eqs. QT%D - ( PH ) and sketched in Fig A 
allows us to calculate the function a' /a. This function is sketched in Fig. 5. In all the 
theoretical generality, the left- hand- side of the barrier in Fig.5 could also consist of several 
pieces, but we do not consider this possibility here. The graph also shows the important 
wave numbers n H , n 2 , n s , ri\. The n H marks the wave whose today's wavelength X(r] R ) = 
2na(r] R )/nH is equal to the today's Hubble radius With our parametrization a(r] R ) = 
21h, this wavenumber is tin = 47r. The n 2 marks the wave whose wavelength A(77 2 ) = 
2na(r]2)/n2 at 77 = 772 is equal to the Hubble radius /(7/ 2 ) at 77 = 772. Since \{j] R ) / \{j]2) = 
(n2/n H )[a(r}R)/a(r)2)] and l{r} R )/l{r} 2 ) = \a(r] R ) / a(i] 2 )}{a{r] R ) / a(r]2)} l/2 , this gives us n 2 /n H = 
[a(?7^)/a(7/ 2 )] 1//2 = (^ 2 - Working out in a similar fashion other ratios, we find 

^ = d, ^ = C, ^ = C>™. (25) 
n H n 2 n s 

Solutions to the gravitational wave equations exist for any a(r]). At intervals of power- 
law dependence a(r]), solutions to Eq. (|TT|) have simple form of the Bessel functions. We 
could have found piece-wise exact solutions to Eq. (|TT|) and join them in the transition 
points. However, we will use a much simpler treatment which is sufficient for our purposes. 
We know that the squeeze parameter r n stays constant in the short-wavelength regimes 
and grows according to Eq. (|TBp in the long-wavelength regime. All modes start in the 
vacuum state, that is, r n = initially. After the end of amplification, the accumulated 
value (|i~7D stays constant up to today. To find today's value of e r " we need to calculate 
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the ratio a**(n) / a*(n) . For every given n, the quantity a* is determined by the condition 
A (r/*) = /(r/*), wheras a** is determined by the condition X(rj^) = l(r]**). 




Let us start from the mode n — n\. For this wave number we have a* = a** = 0(771), and 
therefore r ni = 0. The higher frequency modes n > n\ (above the barrier in Fig.5) have 
never been in the amplifying regime, so we can write 

e r ™ = 1, n > m. (26) 

Let us now consider the modes n in the interval rii > n > n s . For a given n we need to 
know a*(n) and a**(n). Using Eq. flTB) one finds a*(n)/a*(ni) = (ni/n) 1+ ^, and using Eq. 
(|T9|) one finds a^^/a^n,,) = (n s /n) 1+l3s . Therefore, one finds 

a**(n) a**(7i s ) / 77 \ 1+13 



a*{n) a*(77i) V 77 / \77i, 
Since a^{n s ) = a(r] s ), a*(rii) = 0.(771), and a{rj s ) / a{rji) = Ci = ( n i/ n s) 1+l3s , we arrive at 

a** (77) / 77 



a*(77) \77i 

Repeating this analysis for other intervals of the decreasing 77, we come to the conclusion 
that 

e = , 77i > 77 > 77<j, 



77l/ 

e ™ = ( — ) ( — ) > n s >n>n 2 , 



n s j \n x J 
n 



) H M f W^>.>-,, (27) 

77 2 / Vni/ \77i/ 
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The mnemonic rule of constructing e r " at successive intervals of decreasing n is simple. If 
the interval begins at n x , one takes (n/nx) 13 *"^** and multiples with e Tnx , that is, with the 
previous interval's value e Tn calculated at the end of that interval n x . For the function a' /a 
that we are working with, the /3* is always (3, whereas the takes the values (3 S , 0, 1 at 
the successive intervals. 

The modes with n < nu are still in the long-wavelength regime. For these modes, we 
should take a(r]n) instead of ). Combining with a*(n), we find 

e r, = ( jLffafVij'f!:)-*, „<„, (28) 

\n H J \n 2 J Vni/ \n\J 

Formulas (EH) - (|28| ) give approximate values of r n for all n. The factor e Tn is e r ™ > 1 for 
n < ni, and e r ' n ^> 1 for n Cni. This factor determines the mean square amplitude of the 
gravitational waves. 

The mean value of the field hij is zero at every moment of time rj and in every spatial 
point x: (G\hij(r), x)|0) = 0. The variance 

(01^(77, x^'fax^O) = (h 2 ) 

is not zero, and it determines the mean square amplitude of the generated field - the quantity 
of interest for the experiment. Taking the product of two expressions (^[) one can show that 

(h 2 ) = ±-J nJ2h n (v)\dn= h\n, V )™. (29) 
2tx z jo ^ 1 Jo n 

Using the representation (0), (j8|) in Eq. (|29|) one can also write 

(/t 2 ) = / ndn(cosh 2r n + cos 20 n sinh 2r n ). (30) 

We can now consider the present era and use the fact that e Tn are large numbers for all n 
in the interval of our interest n\ > n > nu- Then, we can derive 

h(n,T)) w — / ne r " cos0 w (?7) = 8y/n ( ^ J f— ) e rn cos0„(r]) . (31) 

The quantity h(n,n) is the dimensionless spectral amplitude of the field whose numerical 
value is determined by the calculated squeeze parameter r n . The oscillatory factor cos 4> n {rj) 
reflects the squeezing (standing wave pattern) acquired by modes with n\ > n > n H . For 
modes with n < nn this factor is approximately 1. For high-frequency modes n ^> i%h one 
has 4> n (r]) w n(r] — r] n ) ^> 1, so that h(n,T]) makes many oscillations while the scale factor 
a{rf) is practically fixed at a{rj B ). 

The integral ([30]) extends formally from to 00. Since r n pa for n > n\, the integral 
diverges at the upper limit. This is a typical ultra-violet divergence. It should be discarded 
(renormalized to zero) because it comes from the modes which have always been in their 
vacuum state. At the lower limit, the integral diverges, if /3 < —2. This is an infra-red di- 
vergence which comes from the assumption that the amplification process has started from 
infinitely remote time in the past. One can deal with this divergence either by introduc- 
ing a lower frequency cut-off (equivalent to the finite duration of the amplification) or by 
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considering only the parameters (3 > —2, in which case the integral is convergent at the 
lower limit. It appears that the available observational data (see below) favour this second 
option. The particular case j3 = — 2 corresponds to the de Sitter evolution a{rj) oc \rj\~ 1 . In 
this case, the h(n) found in Eqs. (j3lD, (fiEj) does not depend on n. This is known as the 
Harrison-Zeldvich, or scale- invariant, spectrum. 

An alternative derivation of the spectral amplitude h(n) uses the approximate solutions 
(|T2|), ( |T5D to the wave equation ([□]). This method gives exactly the same, as in Eqs. (^TJ), 
(f26|) - (|28|) numerical values of h(n), but does not reproduce the oscillatory factor cos 4> n (v)- 

One begins with the initial spectral amplitude hi(n) defined by quantum normalization: 
hi(n) = 8y / 7r(Zp;/Aj). This is the amplitude of the mode n at the moment rj* of entering the 
long wavelength regime, i.e. when the mode's wavelength Aj is equal to the Hubble radius 
/(??*). For Aj one derives 

Thus, we have 

hi{n) = A 

where A denotes the constant 



-) , (33) 



A = 68v^F^ . (34) 

The numbers hi(n) are defined at the beginning of the long- wavelength regime. In other 
words, they are given along the left-hand-side slope of the barrier in Fig.5. We want to 
know the final numbers (spectral amplitudes) h(n) which describe the field today, at rjn. 

According to the dominant solution h n (rj) = const of the long- wavelength regime (see 
Eq. (|15D ), the initial amplitude hi(n) stays practically constant up to the end of the long- 
wavelength regime at r/^, that is, up to the right-hand-side slope of the barrier. [The second 
term in Eq. ( |T3"D could be important only at the z-stage and only for parameters (3 S < —(1/2), 
which correspond to the effective equations of state p > e. In order to keep the analysis 
simple, we do not consider those cases.] After the completion of the long-wavelength regime, 
the amplitudes decrease adiabatically in proportion to l/a(r/), up to the present time. Thus, 
we have 

h(n) = A (JL) 2+e ^. (35) 
\n H J a(r] R ) 

Let us start from the lower end of the spectrum, n < uh, and go upward in n. The 
modes n < hh have not started yet the adiabatic decrease of the amplitude, so we have 

/ n \ 2+f3 

h(n) = A — , n<n H . (36) 
\n H J 

Now consider the interval n2 > n > nu- At this interval, the a**(n)/a(?7R) scales as (nn/n) 2 , 
so we have 
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h{n) 



At the interval n s > n > n 2 the ratio a^(n) /a(rjji 
(n 2 /n)(riH/n 2 ) 2 , so we have 



n 2 > n > n H . (37) 
[a**(n) / a(r} 2 )][a(r] 2 ) / a(r} R )] scales as 



M«) = 4-)' + ^, 

\n H J n 2 



n s >n> n 2 . 



Repeating the same analysis for the interval n\ > n > n s we find 

i+P-P* / n \P* n H 



h(n) — A( — 



n 



n H J n 2 



Hi ^ n > n s . 



(38) 



(39) 



It is seen from Eq. ( |3"9"D that an interval of the 2;-stage with (3 S < (the already imposed 
restrictions require also (—1/2) < (3 S ) bends the spectrum h(n) upwards, as compared with 
Eq. (PS]), for larger n. If one recalls the relationship (^) between l Q and Ih and uses 
(p8|) in Eq. ([IT]) one arrives exactly at Eqs. (|36|)-(|39D up to the oscillating factor cos (j) n (rj). 

Different parts of the barrier in Fig.5 are responsible for amplitudes and spectral slopes 
at different intervals of n. The sketch of the generated spectrum h(n) in conjunction with 
the form of the barrier is shown in Fig. 6. 




Fig. 6. Amplitudes and spectral slopes of h(n) are determined by different 

parts of the barrier a '/a. 



The present day frequency of the oscillating modes, measured in Hz, is defined as v = 
cn/2na(r]R). The lowest frequency (Hubble frequency) is vu = c/Ih- For numerical estimates 
we will be using vu w 10~ 18 Hz. The ratios of n are equal to the ratios of z/, so that, for 
example, n/n^ = vjvu. For high-frequency modes we will now often use the ratios of v 
instead of ratios of n. 
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In addition to the spectral amplitudes h{n) the generated field can be also characterized 
by the spectral energy density parameter Q g (n). The energy density e g of the gravitational 
wave field is 

Kt 9 = ~^i\ h ijfi = -^h tj hi/. 
The mean value (0|e ff (r/, x)|0) is given by 



1 C 2 f°° 2 \s> 



o 



s=l 



For high-frequency modes, it is only the factor e ±inr> that needs to be differentiated by 77. 

s 

After avaraging out the oscillating factors, one gets h 



n 2 



fh 



so that 



^ = if" 2ft2 w^ < 41 > 

In fact, the high-frequency approximation, that has been used, permits integration over 
lower n only up to uh- And the upper limit, as was discussed above, is in practice n\, not 
infinity The parameter Q g is defined as Q g = (e g )/e, where e is given by Eq. ( |2"3"D (critical 
density). So, we derive 

rni r[ n rv\ Ay 

n g = / n g (n)— = / n g (p)— 

J Tin J uh V 

and 



V ^ 2 



%{v) = -h\v)[—J . (42) 

The dimensionless quantity Vt g {u) is useful because it allows us to quickly evaluate the 
cosmological importance of the generated field in a given frequency interval. However, 
the primary and more universal concept is h(v), not Vl g {v). It is the field, not its energy 
density, that is directly measured by the gravity-wave detector. One should also note that 
some authors use quite a misleading definition fl g (f) = (1/ p c )(dp gw / din f) which suggests 
differentiation of the gravity-wave energy density by frequency. This would be incorrect and 
could cause disagreements in numerical values of Q g . Whenever we use fi 9 (z/), we mean 
relationship (|42|); and for order of magnitude estimates one can use |L|: 

n g (u)^h 2 (u)(—) 2 . (43) 



V. THEORETICAL AND OBSERVATIONAL CONSTRAINTS 

The entire theoretical approach is based on the assumption that a weak quantized gravity- 
wave field interacts with a classical pump field. We should follow the validity of this approx- 
imation throughout the analysis. The pump field can be treated as a classical gravitational 
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field as long as the driving energy density e is smaller than the Planck energy density, or, in 
other words, as long as the Hubble radius l(rj) is greater than the Planck length Ipi. This 
is a restriction on the pump field, but it can be used as a restriction on the wavelength Aj 
of the gravity- wave mode n at the time of entry the long- wavelength regime. If /(r/*) > Ipi, 
then Aj > Ipi. The A, is given by Eq. fl32|). So, we need to ensure that 



lo \VhJ 

At the lowest-frequency end v — i/jj this inequality gives b(lpi/l ) < 1. In fact, the observa- 
tional constraints (see below) give a stronger restiction: 

b ki ~ 10" 6 , (44) 
which we accept. Then, at the highest-frequency end v = v\ we need to satisfy 

< 10 6 . (45) 



Let us now turn to the generated spectral amplitudes h(v). According to Eq. fl36D 
we have h{i>p) « 68y / 7r(/pi// ). The measured microwave beckgound anisotropics, which 
we discuss below, require this number to be at the level of 10~ 5 , which gives the already 
mentioned Eq. (f4l|). The quantity h{v\) at the highest frequency V\ is given by Eq. ([39]): 

Using Eq. fl22|) this expression for hiv\) can be rewritten as 

I'M Vp Al 

where Ai = c/u%. This last expression for h{v{) is not surprising: the modes with v > V\ 
are still in the vacuum state, so the numerical value of h{i>\) is determined by quantum 
normalization. 

All the amplified modes have started with small initial amplitudes hi, at the level of 
zero-point quantum fluctuations. These amplitudes are also small today, since the hi could 
only stay constant or decrease. However, even these relatively small amplitudes should obey 
observational constraints. We do not want the fl g in the high-frequency modes, which might 
affect the rate of the primordial nucleosynthesis, to exceed the level of 10 -5 . This means 
that Q g (i/i) cannot exceed the level of 10~ 6 or so. The use of Eq. (E2j) in combination 
with Qg(ui) w 10~ 6 and h{i>\) from Eq. (f46|) , gives us the highest allowed frequency v\ ~ 
3 x 10 10 Hz. We will use this value of v\ in our numerical estimates. Returning with this 
value of V\ to Eq. (|45| ) we find that parameter {3 can only be (3 < —1.8. We will be treating 
(3 = —1.8 as the upper limit for the allowed values of (3. 

We can now check whether the accepted parameters leave room for the postulated z-stage 
with (3 S < 0. Using Eq. ( p2|) we can rewrite Eq. (|44] ) in the form 
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We know that v^jvn — 10 2 and V\jv s is not smaller than 1. Substituting all the numbers in 
Eq. ([47D one can find that this equation cannot be satisfied for the largest possible (3 = —1.8. 
In the case f3 = —1.9, Eq. ( |4"7| ) is only marginally satisfied, in the sense that a significant 
deviation from (3 S = toward negative (3 S can only last for a relatively short time. For 
instance, one can accomodate (3 S = —0.4 and v s = 10 s Hz. On the other hand, if one takes 
(3 = —2, a somewhat longer interval of the z-stage with (3 S < can be included. For instance, 
Eq. (|47|) is satisfied if one accepts v s = 10~ 4 Hz and (3 S = —0.3. This allows us to slightly 
increase h(i>) in the interval v s < v < b>i, as compared with the values of hiy) reached in 
the more traditional case (3 = —2, f3 s = 0. In what follows, we will consider consequences 
of this assumption for the prospects of detection of the produced gravitational wave signal. 
Finally, let us see what the available information on the microwave background anisotropies 
|T5| |, [|nj allows us to conclude about the parameters /3 and l a . 

Usually, cosmologists operate with the spectral index n (not to be confused with the wave 
number n) of primordial cosmological perturbations. Taking into account the way in which 
the spectral index n is defined, one can relate n with the spectral index (3 + 2 that shows up 
in Eq. ([36]). The relationship between them is n = 2(3 + 5. This relationship is valid inde- 
pendently of the nature of cosmological perturbations. In particular, it is valid for density 
perturbations, in which case the h(n) of Eq. ( j36|) is the dimensionless spectral amplitude 
of metric perturbations associated with density perturbations. If primordial gravitational 
waves and density perturbations were generated by the mechanism that we discuss here (the 
assumption that is likely to be true) than the parameter (3 that participates in the spectral 
index is the same one that participates in the scale factor of Eq. (0). Primordial gravi- 
tational waves and primordial density perturbations with the same spectral index produce 
approximately the same lower-order multipole distributions of large-scale anisotropies. 

The evaluation of the spectral index n of primordial perturbations have resulted in n = 
1.2 ±0.3 [|1(J or even in a somewhat higher value. A recent analysis [H] of all available data 
favors n = 1.2 and the quadrupole contribution of gravitational waves twice as large as that 
of density perturbations. One can interpret these evaluations as indication that the true 
value of n lies somewhere near n = 1.2 (hopefully, the planned new observational missions 
will determine this index more accurately). This gives us the parameter (3 somewhere near 
(3 = —1.9. We will be using (3 = —1.9 in our estimates below, as the observationally preferred 
value. The parameter (3 can be somewhat larger than (3 = —1.9. However, as we already 
discussed, the value (3 = —1.8 (n = 1.4) is the largest one for which the entire approch is 
well posed. The Harrison-Zeldovich spectral index n = 1 corresponds to (3 = —2. 

The observed quadrupole anisotropy of the microwave background radiation is at the 
level 5T/T « 10~ 5 . The quadrupole anisotropy that would be produced by the spectrum 
(f36|) - ([39]) is mainly accounted for by the wave numbers near n#. Thus, the numerical value 
of the quadrupole anisotropy produced by relic gravitational waves is approximately equal 
to A. According to general physical considerations and detailed calculations |18[], the metric 
amplitudes of long-wavelength gravitational waves and density perturbations generated by 
the discussed amplification mechanism are of the same order of magnitude. Therefore, they 
contribute roughly equally to the anisotropy at lower multipoles. This gives us the estimate 
A fa 10~ 5 , that we have already used in Eq. (44). It is not yet proven observationally that a 
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significant part of the observed anosotropies at lower multipoles is indeed provided by relic 
gravitational waves, but we can at least assume this with some degree of confidence. It is 
likely that the future measurements of the microwave background radiation will help us to 
verify this theoretical conclusion. 

Combining all the evaluated parameters together, we show in Fig.7 the expected spec- 
trum of h(y) for the case f3 = —1.9. A small allowed interval of the z-stage is also included. 
The intervals of the spectrum accessible to space-based and ground-based interferometers 
are indicated by vertical lines. 



log h(v) 




logv, Hz 

Fig. 7. Expected spectrum h(v) for the case (3 = -1.9. 



It is necessary to note [18|], |19[ that the confirmation of any n > 1 (/3 > —2) would mean 



that the very early Universe was not driven by a scalar field - the cornerstone of inflationary 
considerations. This is because the n > 1 (/3 > —2) requires the effective equation of state at 
the initial stage of expansion to be e + p < (see Eq. (H)), but this cannot be accomodated 
by any scalar field with whichever scalar field potential. The available data do not prove yet 
that n > 1, but this possibility seems likely. 

It is also necessary to say that a certain damage to gravitational wave research was 
inflicted by the so called "standard inflationary result" . The "standard inflationary result" 
predicts infinitely large amplitudes of density perturbations in the interval of spectrum with 
the Harrison- Zeldovich slope n = 1 (J3 — —2): Sp/p oc 1 fy/l — n. The metric (gravitational 
field) amplitudes of density perturbations are also predicted to be infinitely large, in the 
same proportion. Through the so-called "consistency relation" this divergence leads to the 
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vanishingly small amplitudes of relic gravitational waves. Thus, the "standard" inflationary 
theory predicts zero for relic gravitational waves; the spectrum similar in shape to the one 
shown in Fig. ,7 would have been shifted down by many orders of magnitude. This prediction 
is hanging on the "standard inflationary result" , but the "result" itself is in a severe conflict 
not only with theory but with observations too: when the observers marginalize their data to 
n = 1 (enforce this value of n in data analysis) they find finite and small density perturbations 
instead of infinitely large perturbations predicted by inflationary theorists. [For analytical 
expressions of the "standard inflationary result" see any inflationary article, including recent 
reviews. For graphical illustration of the predicted divergent density perturbations and 
quadrupole anisotropies see, for example, [p0| . For critical analysis and disagreement with 
the "standard inflationary result" see ||18|| .1 General relativity and quantum field theory do 
not produce the "standard inflationary result" , so we shall better return to what they say. 



VI. DETECTABILITY OF RELIC GRAVITATIONAL WAVES 

We switch now from cosmology to prospects of detecting the predicted relic gravitational 
waves. The ground-based PTJ- P3[ and space-based El, [EH] laser interferometers (see also 



|26| - |p8[| ) will be in the focus of our attention. We use laboratory frequencies v and intervals 
of laboratory time t (cdt = airjRjdrf). Formulas fl38|) and (|39|) , with A = 1CT 5 , z^/Vh = 10 2 , 
and the oscillating factor restored, can be written as 



h{v, t) « 1(T 7 cos[27rz/(t - i„)] — , v 2 < v < v s (48) 



and 



( v \ i+z 3 -/^ f v \@ s 

h{v, t) « KT 7 cos[2vrz/(t - £„)] — — . u, < v < v x (49) 

\v H ) \Vh' 

where the deterministic (not random) constant t u does not vary significantly from one fre- 
quency to another at the intervals Au « v. The explicit time dependence of the spectral 
variance h 2 {v, t) of the field, or, in other words, the explicit time dependence of the (zero-lag) 
temporal correlation function of the field at every given frequency, demonstrates that we are 
dealing with a non-stationary process (a consequence of squeezing and severe reduction of 
the phase uncertainty). We will first ignore the oscillating factor and will compare the pre- 
dicted amplitudes with the sensitivity curves of advanced detectors. The potential reserve 
of improving the signal to noise ratio by expoloiting the squeezing will be discussed later. 

Let us start from the Laser Interferometer Space Antenna (LISA) [p4|l . The instrument 
will be most sensitive in the interval, roughly, from 10~ 3 Hz to 10 _1 ifz, and will be rea- 
sonably sensitive in a broader range, up to frequencies 10~ A Hz and 1Hz. The sensitivity 
graph of LISA to a stochastic background is usually plotted under the assumption of a 
1-year observation time, that is, the root-mean-square (r.m.s.) instrumental noise is being 
evaluated in frequency bins Au = 3 x 10~ 8 Hz around each frequency v. We need to rescale 
our predicted amplitude h{y) to these bins. 

The mean square amplitude of the gravitational wave field is given by the integral (PU[). 
Thus, the r.m.s. amplitude in the band Av centered at a given frequency v is given by the 
expression 
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h(p,Av) = h(v)^J—. (50) 

We use Eqs. (|3J), ( f49"D and calculate expression fl5"0|) assuming Au = 3xlO~ 8 Hz. The results 
are plotted in Fig.8. Formula fl48|) has been used throughout the covered frequency interval 
for the realistic case (3 = —1.9 and for the extreme case (3 = —1.8. The line marked z-model 
describes the signal produced in the composite model with (3 = — 2 up to v s = 10~ 4 Hz 
(formula (pE8"D) and then followed by formula ( pE9"D with (3 S = —0.3. This model gives the 
signal a factor of 3 higher at v — 10~~ 3 Hz, than the model f3 = —2 extrapolated down to 
this frequency. 



LISA Sensitivity 




10" 4 10" 3 10" 2 10" 1 10° 

frequency V (Hz) 



Fig. 8. Expected spectrum (3 = -1.9 and other possible spectra 
in comparison with the LISA sensitivity. 



There is no doubt that the signal (3 = —1.8 would be easily detectable even with a single 
instrument. The signal (3 = —1.9 is marginally detectable, with the signal to noise ratio 
around 3 or so, in a quite narrow frequency interval near and above the frequency 3x 10~ 3 Hz. 
However, at lower frequencies one would need to be concerned with the possible gravitational 
wave noise from unresolved binary stars in our Galaxy. The further improvement of the 
expected LISA sensitivity by a factor of 3 may prove to be crucial for a confident detection 
of the predicted signal with (3 = —1.9. 

Let us now turn to the ground-based interferometers operating in the interval from 10Hz 
to 10 4 Hz. The best sensitivity is reached in the band around v = 10 2 Hz. We take this 
frequency as the representative frequency for comparison with the predicted signal. We will 
work directly in terms of the dimensionless quantity hiy). If necessary, the r.m.s. amplitude 
per Hz 1 / 2 at a given v can be found simply as h{y)j^fv. The instrumental noise will also 
be quoted in terms of the dimensionless quantity h ex {y). 
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The expected sensitivity of the initial instruments at v — 10 2 Hz is h ex = 10 -21 or 
better. The theoretical prediction at this frequency, following from (fi8|), d49|) with f3 s = 0, 
is h t h = 1CT 23 for (3 = —1.8, and h t h = 10~ 25 for /3 = —1.9. Therefore, the gap between 
the signal and noise levels is from 2 to 4 orders of magnitude. The expected sensitivity of 
the advanced interferometers, such as LIGO-II fl29j, can be as high as h ex = 1CT 23 . In this 



case, the gap vanishes for the (3 = —1.8 signal and reduces to 2 orders of magnitude for 
the j3 = —1.9 signal. Fig.9 illustrates the expected signal in comparison with the LIGO-II 
sensitivity. Since the signal lines are plotted in terms of h(v), the LISA sensitivity curve 
(shown for periodic sources) should be raised and adjusted in accordance with Fig.S. 




A signal below noise can be detected if the outputs of two or more detectors can be 
cross correlated. [For the early esimates of detectability of relic gravitational waves see [50].] 
The cross correlation will be possible for ground-based interferometers, several of which are 
currently under construction. The gap between the signal and the noise levels should be 
covered by a sufficiently long observation time r. The duration r depends on whether the 
signal has any temporal signature known in advance, or not. We start from the assumption 
that no temporal signatures are known in advance. In other words, we first ignore the 
squeezed nature of the relic background and work under the assumption that the squeezing 
cannot be exploited to our advantage. 

The response of an instrument to the incoming radiation is s(t) = F^h 1 ^ where F^ 
depends on the position and orientation of the instrument. Since the is a quantum- 
mechanical operator (see Eq. (Bf)) we need to calculate the mean value of a quadratic 
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quantity. The mean value of the cross correlation of responses from two instruments 
(0|si(t)s 2 (£)|0) will involve the overlap reduction function |3lJ|- [Q, which we assume to 
be not much smaller than 1 ||33|| . The signal to noise ratio S/N in the measurement of the 
amplitude of a signal with no specific known features increases as (rz/) 1 / 4 , where v is some 
characteristic central frequency. If the signal has features known in advance and exploited 
by the matched filtering technique, the S/N increases as (rz/) 1 / 2 . 

We apply the guaranteed law (rz/) 1 / 4 to initial and advanced instruments at the repre- 
sentative frequency v = 10 2 Hz. This law requires a reasonably short time r = 10 6 sec in 
order to improve the S/N in initial instruments by two orders of magnitude and to reach 
the level of the signal with extreme spectral index j3 = —1.8. The longer integration time or 
a better sensitivity will make the S/N larger than 1. In the case of a realistic spectral index 
(3 = —1.9 the remaining gap of 4 orders of magnitude can be covered by the combination 
of a significantly better sensitivity and a longer observation time (not necessarily in one 
non- interrupted run). The sensitivity of the advanced laser interferometers, such as LIGO 
II, at the level h ex = 10 -23 and the same observation time r = 10 6 sec would be sufficient 
for reaching the level of the predicted signal with /3 = —1.9. 

An additional increase of S/N can be achieved if the statistical properties of the signal 
can be properly exploited. Squeezing is automatically present at all frequencies from v H 
to Ui. The squeeze parameter r is larger in gravitational waves of cosmological scales, and 
possibly the periodic structure in Eq. ([H]) can be better revealed at those scales. However, 
we are interested here in frequencies accessible to ground based interferometers, say, in the 
interval 30Hz — 100Hz. If our intention were to monitor one given frequency v from the 
beginning of its oscillating regime and up till now, then, in order to avoid the destructive 
interference from neighbouring modes during all that time, the frequency resolution of the 
instrument should have been increadibly narrow, of the order of 10~ 18 Hz. Certainly, this 
is not something what we can, or intend to do. Although the amplitudes of the waves 
have adiabatically decreased and their frequencies redshifted since the beginning of their 
oscillating regime, the general statistical properties of the discussed signal are essentially 
the same now as they were 10 years after the Big Bang or will be 1 million years from now. 

The periodic structure (^) may survive at some level in the instrumental window of 
sensitivity from v min (minimal frequency) to v max (maximal frequency). The mean square 
value of the field in this window is 

/ h i^ t )— = 10- u ^— cos 2 [27rz/(t-t,)]^ +1 dz/. (51) 

Because of the strong dependence of the integrand on frequency, v~ 2S or z/~ 2 ' 8 , the value of 
the integral ([51]) is determined by its lower limit. Apparently, the search through the data 
should be based on the periodic structure that may survive at v = v min . As an illustration, 
one can consider such a narrow interval Av = Umax — ^min that the integral ([)!]) can be 
approximated by the formula 

/ h (u, t) — w 10 cos [2iris min (t - t min )\ . 

J Vmin 

Clearly, the correlation function is strictly periodic and its structure is known in advance, in 
contrast to other possible signals. This is a typical example of using the appriori information. 
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Ideally, the gain in S/N can grow as {Tv min ) l l 2 . This would significantly reduce the required 
observation time r. For a larger Az/, even an intermediate gain between the guaranteed 
law (rz/) 1//4 and the law (rz/) 1 / 2 , adequate for the matched filtering technique, would help. 
This could potentially make the signal with f3 = — 1.9 measurable even by the initial laser 
interferometers. A straightforward application of ( pi]) for exploiting the squeezing may not 
be possible, as argued in the recent study p5[ , but more sophisticated methods are not 
excluded. 

For frequency intervals covered by bar detectors and electromagnetic detectors, the ex- 
pected results follow from the same formulas (08|), (f49[) and have been briefly discussed 



elsewhere 30 , IE 



VII. CONCLUSION 



It would be strange, if the predicted signal at the level corresponding to (3 — —1.9 were 
not seen by the instruments capable of its detection. There is not so many cosmological 
assumptions involved in the derivation, that could prove wrong, thus invalidating our pre- 
dictions. On the other hand, it would be even more strange (and even more interesting) if 
the relic gravitational waves were detected at the level above the (3 = — 1.8 line. This would 
mean that there is something fundamentally wrong in our basic cosmological premises. To 
summarise, it is quite possible that the detection of relic (squeezed) gravitational waves 
may be awaiting only the first generation of sensitive instruments and an appropriate data 
processing strategy. 



VIII. ACKNOWLEGEMENTS 



I appreciate the help of M. V. Prokhorov in preparation of the figures. 



[1] L. P. Grishchuk, Zh. Eksp. Teor. Fiz. 67, 825 (1974) [JETP 40, 409 (1975)]; Ann. NY Acad. 

Sci. 302, 439 (1977); Uspekhi Fiz. Nauk. 156, 297 (1988) [Sov. Phys.-Uspekhi. 31, 940 (1988)]. 
[2] K. S. Thorne, in 300 Years of Gravitation, eds. S. W. Hawking and W. Israel (Cambridge: 

CUP) 1987, p. 330. 

[3] K. S. Thorne, in Particle and Nuclear Astrophysics and Cosmology in the Next Millenium, 

eds. E. Kolb and R. Peccei (Singapore: World Scientific) 1995, p. 160. 
[4] B. F. Schutz, Class. Quant. Grav. 16, A131 (1999). 

[5] P. L. Knight, in Quantum Fluctuations, Eds. S. Reynaud, E. Giacobino, and J. Zinn-Justin, 

(Elsevier Science) 1997, p. 5. 
[6] L. P. Grishchuk and Yu. V. Sidorov, Class. Quant. Grav. 6, L161 (1989); Phys. Rev. D42, 

3413 (1990). 

[7] L. P. Grishchuk, in Workshop on Squeezed States and Uncertainty Relations, NASA Conf. 
Publ. 3135, 1992, p. 329; Class. Quant. Grav. 10, 2449 (1993); in Quantum Fluctuations, 
Eds. S. Reynaud, E. Giacobino, and J. Zinn-Justin, (Elsevier Science) 1997, p. 541. 

[8] L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields (New York: Pergamon) 1975. 



24 



W. Schleich and J. A. Wheeler, J. Opt. Soc. Am B4, 1715 (1987); W. Schleich et. al, Phys. 
Rev. A40, 7405 (1989). 

Ya. B. Zeldovich and I. D. Novikov, The Structure and Evolution of the Universe (Chicago, 
IL: University of Chicago Press) 1983. 

M. Giovannini, Phys. Rev. D58, 083504 (1998); Report |hep^ph/991248C . 
L. P. Grishchuk, Report gr-qc/9810055[ 

G. Veneziano, Phys. Lett. B265, 287 (1991): M. Gasperini and G. Veneziano, Astropart. Phys. 

1, 317 (1993); M. Gasperini and M. Giovannini, Phys. Rev. D47, 1519 (1993); M. Gasperini, 

Report |hep-th/9607146 . 

T. Creighton, Report gr-qc/9907045 . 

G. F. Smoot et. al, Astroph. J. 396, LI (1992). 

C. L. Bennet et. al, Astroph. J. 464, LI (1996). 

A. Melchiori, M. V. Sazhin, V. V. Shulga, and N. Vittorio, Astroph. J. 518, 562 (1999). 
L. P. Grishchuk, Phys. Rev. D50, 7154 (1994); in Current Topics in Astro fundamental Physics: 
Primordial Cosmology, Eds. N. Sanchez and A. Zichichi, (Kluwer Academic) 1998, p. 539; 
Report |gr-qc/98010TT. 



L. P. Grishchuk, Class. Quant. Grav. 14, 1445 (1997). 
J. Martin and D. J. Schwarz, Report astro-ph /991 1225| . 

A. Abramovici et. al, Science 256, 325 (1992). 

C. Bradaschia et. al, Nucl. Instrum. and Methods A289, 518 (1990). 
J. Hough and K. Danzmann et. al, GEO600 Proposal, 1994. 
P. Bender et. al. LISA Pre-Phase A Report, Second Edition, 1998 
S. L. Larson, W. A. Hiscock. and R. W. Hellings, Report |gr^qc/990908C . 
Gravitational Wave Experiments, eds. E. Coccia, G. Pizzella, and F. Ronga (Singapore: World 
Scientific) 1995. 

Gravitational Waves: Sources and detection, eds. I. Giufolini and F. Fidecaro (Singapore: 
World Scientific) 1997. 

Laser Interferometer Space Antenna, ed. W. Folkner (AIP Conference Proceedings 456) 1998. 
LIGO II Conceptual Project Book, LIGO-M990288-00-M. 

L. P. Grishchuk, Pis'ma Zh. Eks. Teor. Fiz. 23, 326 (1976) [JETP Lett. 23, 293 (1976)]. 
P. F. Michelson, Mon. Not. R. Astr. Soc. 227, 933 (1987). 
N. L. Christensen, Phys. Rev. D46, 5250 (1992). 
E. E. Flanagan, Phys. Rev. D48, 2389 (1993). 

B. Allen, in Relativistic Gravitation and Gravitational Radiation, Eds. J-A. Marck and J-P. 
Lasota (CUP, 1997) p. 373. 

[35] B. Allen, E. E. Flanagan, and M. A. Papa. Report [gTqc /9906054 . 



25 



